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The influence of coherent optical nonlinearities on polariton propagation effects is studied within 
a theory-experiment comparison. A novel approach that combines a microscopic treatment of the 
boundary problem in a sample of finite thickness with excitonic and biexcitonic nonlinearities is 
introduced. Light-polarization dependent spectral changes are analyzed for single-pulse transmission 
and pump-probe excitation. 

PACS numbers: 71.36.+c,71.35.Cc 
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The propagation of light pulses which are coupled to 
the excitonic resonances is a fundamental problem in 
semiconductor optics that has been the subject of intense 
experimental and theoretical research. The observations 
are dominated by the strong light-matter interaction and 
its interplay with the inherent many-particle Coulomb in- 
teraction of the electronic system. In the linear optical 
regime, exciton polaritons give rise to effects like polari- 
ton beating in the time-resolved pulse transmission [lj or 
yield strong modifications of excitonic transmission spec- 
tra, see e.g. In the nonlinear optical regime, in- 
coherent saturation of the polariton resonances in trans- 
mission spectra 0, and in the amplitude and phase 
of the time-resolved transmission as well as transmis- 
sion changes in pump-probe experiments have been 
studied. 

The above investigations are focused on samples where 
the thickness corresponds to a few exciton Bohr radii 
where polariton effects are most pronounced in the op- 
tical spectra. In this case the interplay of the induced 
excitonic polarization in the medium with the propagat- 
ing light field is strongly influenced by sample surfaces 
which considerably complicates the theoretical descri 
tion. In a frequently used phcnomcnological approach 
the coupling of the exciton relative and center-of-mass 
(COM) motion at the sample boundaries is neglected. 
Then the exciton COM motion is subject to quantiza- 
tion effects in the confinement geometry while the ex- 
citon relative motion is approximated by the result of 
the infinitely extended medium. As a result of this ap- 
proach, the solution of the wave equation for the electro- 
magnetic field requires additional boundary conditions 
(abcs) HE3 which are, however, not uniquely defined. 
Unfortunately, in many situations the choice of ABCs 
can influence the theoretical predictions To avoid 

these ambiguities, the use of microscopic boundary con- 



ditions within a non-local semiconductor response has 
been discussed in Refs. |l2 l lF j| and recently applied to 
the linear optical regime p&ETIlil. The description of 
optical nonlinearities combined with a microscopic treat- 
ment of boundary conditions has been restricted so far 
to the quantum-well (QW) limit where propagation ef- 
fects lead to radiative exciton broadening, which can be 
modified in radiatively coupled QWs. Furthermore, ex- 
citonic nonlinearities in QWs have been studied in con- 
nection with microcavity polaritons; for a review see |l5j . 
In the opposite limit of the sample thickness approac hing 
the bulk limit nonlinear pulse propagation effects [la, Il7| 
have been analyzed successfully within a local approach 
for the semiconductor response. The decoupling of rela- 
tive and COM exciton motion has also been used in an 
earlier approach to study the influence of propagation 
effects on four- wave mixing signals [lj| • 

In this letter we present a novel approach that de- 
scribes the combined influence of polariton propagation 
effects and coherent excitonic and biexcitonic nonlineari- 
ties. While the theoretical understanding of the coherent 
nonlinear excitation dynamics in QW or bulk systems is 
well developed, the additional influence of propagation 
effects is missing especially in a regime where the sample 
thickness exceeds the QW limit but the heterostructure 
interfaces prevent bulk-like behavior. In this case typi- 
cally several well-resolved polariton modes contribute to 
the optical properties and, as a result of our investiga- 
tions, the Coulomb interaction of the corresponding ex- 
citonic states modifies optical nonlinearities. 

Our approach is based on a direct solution of the two- 
and four-particle Schrodinger equations for the exciton 
and biexciton motion together with Maxwell's equations. 
We apply microscopic boundary conditions to avoid am- 
biguities due to ABCs in situations where both material 
polarization as well as optical fields are strongly influ- 
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enced by the boundaries of the system. A direct compar- 
ison of nonlinear transmission spectra for a ZnSe/ZnSSe 
heterostructure shows excellent agreement between the- 
ory and experiment and supports the detailed analysis 
of the microscopic model. As another application, cal- 
culated pump-probe spectra for a GaAs/AlGaAs het- 
erostructure are presented. 

In the coherent regime, excitonic and biexcitonic non- 
linearities up to third order in the optical field can 
be consistently described in terms of the dynamics- 
controlled truncation (DCT) formalism [l^ . Eof which 
has been successfully applied in the past to QW systems 
[H iHH 0111 El- We extend this formalism to lay- 
ers with a finite thickness where the sample boundaries 
still provide a confinement potential for the electrons and 
holes. To account for the coupling of the exciton relative 
and COM motion, independent coordinates for the elec- 
tron and hole motion in propagation direction are used 
while these carriers are subject to Coulomb interaction 
(treated on the level of DCT) and to the optical field. 

We consider a two-band model with spin-degenerate 
conduction and valence bands to describe the spectrally 
lowest interband transitions for a semiconductor layer in 
a slab geometry jjjj with strain-split light- and heavy- 
hole valence bands. The resonant contribution to the 
spatially inhomogeneous polarization in the semiconduc- 
tor, 
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is given in terms of the non-local excitonic transition 
amplitude pf^ z Zb \(t) — (ft-k(zh)ek(ze)) with the elec- 
tron ek(ze) and hole /ik(^h) operators, respectively. Ac- 
cording to the dipole selection rules, the (— i,— |) and 
transitions are independently driven by circu- 
lar optical polarization with unit vectors e + and e_, 
respectively, k is the in-plane carrier momentum (for 
the directions possessing translational invariance) and 
e, h are quantum numbers that simultaneously denote 
the band index and the z-component of the correspond- 
ing total angular momenta. A local dipole-interaction 
d Gh (k, z c - z h ) = d ch S(z c - z h ) is used. Applying the 
DCT scheme, an equation of motion for the excitonic 
transition amplitude p?£ z > is obtained which is cou- 
pled to the biexcitonic correlation function, 

-(h' kl (ziK 2 (z 2 ))(h k3 (z 3 )e ki (z 4 )) 

+('4 1 (2i)e k4 (*4))(fcka (*3K 2 (z 2 )> , (2) 

that itself obeys a four-particle Schrodinger equation 
with coherent source terms 27] . The resulting closed 
set of equations contains all mean-field (Hartree-Fock) 
as well as biexcitonic contributions to the semiconductor 
response up to third order in the optical field. The eval- 
uation of these material equations is performed in the 



exciton eigenbasis. As an important ingredient of our 
approach, exciton eigenfunctions </> OT (k, z e , Zh) are deter- 
mined for the given confinement geometry to include the 
additional non-local space dependence. In contrast to 
Ref. 01 we use eigenfunctions that individually fulfill the 
boundary conditions of the system 0, [2^ . Then every 
eigenfunction corresponds to a polariton resonance in the 
linear optical spectrum and a truncation of the expansion 
can be justified with the considered spectral range while 
the boundary conditions of the inhomogeneous system re- 
main fully satisfied. The excitonic transition amplitude 
and the singlet (A = —1) and triplet (A = +1) contribu- 
tions to the biexcitonic correlation function are expressed 
as 
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where p^{t) and b^ hx (q,t) are the time-dependent 
expansion coefficients, a = m£/M* and f3 = m* e /M* 
are the ratio of the hole and the electron masses to the 
total exciton mass M* = m* + m£, respectively. This 
treatment generalizes the eigenfunction expansion in the 
QW limit discussed in |2^| . For a self-consistent descrip- 
tion of the light propagation in a semiconductor material, 
Maxwell's equations are solved directly together with 
the excitonic and the biexcitonic dynamics. 

In Fig. H the theoretical and experimental results are 
shown for the transmission of a single light pulse through 
a high-quality 20 nm ZnSe layer sandwiched between two 
ZnSSe cladding layers. The sample was pseudomorphi- 
cally grown by molecular-beam epitaxy on a (OOl)GaAs 
substrate. Further details on growth and sample char- 
acterization as well as on modeling of the finite-height 
confinement potentials at the ZnSe/ZnSSe interfaces are 
described in Q. The opaque GaAs substrate was care- 
fully removed by grinding and subsequent chemical etch- 
ing to permit measurements in transmission geometry. 
The sample was mounted onto a glass plate and kept on 
a cold finger in a microcryostat at a temperature of 4 K. 
A frequency-doubled self-mode-locked Ti:sapphire laser 
with 110 fs (FWHM) pulse duration and 83 MHz repeti- 
tion rate was used as excitation source for the spectrally 
resolved but time integrated transmission experiments. 
The maximum of the spectral pulse profile was set close 
to the spectral position of the first heavy-hole (hhl) po- 
lariton mode to excite the whole spectral region from the 
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FIG. 1: (a) Dashed-dotted line: Linear transmission T (de- 
picted as 1 — T) for a pulse energy of 1.2 pj for the 20 nm 
ZnSe sample. Dotted line: Spectral shape of the 110 fs laser 
pulse. Nonlinear transmission for a pulse energy of 12.2 pj 
for linear (solid line) and circular (dashed line) light polar- 
ization, (b) Optically induced transmission changes AT for 
linear (solid line) and circular (dashed line) light polariza- 
tion, corresponding to the data in (a), (c), (d) Theoretical 
results corresponding to (a) and (b), without the influence of 
Fabry-Perot effects (see discussion in the text). 

exciton-bicxciton transition up to the hh2 and hh3 po- 
lariton resonances as shown by the dotted line in Fig.^i. 
Linear or circular light polarization was selected using a 
Pockels cell. The transmitted signal was analyzed with 
a spectrometer and was recorded with a CCD camera. 

The measured linear transmission spectrum in the 
vicinity of the heavy-hole (hh) exciton resonance, de- 
picted as dashed-dotted line in Fig. exhibits three 
pronounced polariton resonances denoted by hhl-hh3. 
With the microscopic treatment of boundary conditions, 
these polariton modes are perfectly reproduced in Fig. 
(without adjustable parameters like a dead-layer thick- 
ness, the sample thickness has been independently deter- 
mined using x-ray diffraction Q). Results for nonlinear 
transmission spectra are shown in Figs.^i (experiment) 
and c (theory) where linear and circular light polarization 
has been used. For better comparison, the optically in- 
duced transmission changes AT = T non ii near — Tii near are 
given in Figs. ^ an d d, for the experiment and theory, 
respectively. The outermost surfaces between cladding 
layer and air lead to a decreased transmission as well as to 



weak Fabry-Perot effects in the experiment |4| . In order 
to avoid additional parameters, the theory only includes 
intrinsic polariton effects. The constant offset in the mea- 
sured nonlinear transmission results from slight intensity 
drifts of the applied laser pulses. These effects are of mi- 
nor importance for the discussion here, since they just 
cause a nearly constant offset in the investigated trans- 
mission changes. The appearance of a bound biexciton 
(XX) resonance for linear optical polarization, which is 
absent for circular optical polarization according to its 
spin-singlet symmetry, is clearly identified in the theory- 
experiment comparison. The spectral changes at the po- 
lariton resonances result from mean-field (Hartree-Fock) 
as well as from biexcitonic correlations. For circular light 
polarization only the biexcitonic continuum states with 
electronic triplet configuration are excited whereas for 
linear light polarization both singlet and triplet biexci- 
tonic states contribute to the transmitted signal. In par- 
ticular, an advantage of the theoretical approach applied 
here is the simultaneous inclusion of the bound biexciton 
state and the exciton-exciton scattering continuum. The 
latter one is essential to reproduce the broad background 
on the high-energy side of the polariton resonances in the 
nonlinear transmission spectra. As a minor aspect, the 
biexciton binding energy is slightly underestimated by 
the truncation of the exciton basis, similar to the result 
in |26( for a QW system. 

To provide further insight into the nature of coherent 
nonlinear polariton saturation, we also studied the pump- 
probe excitation for a GaAs layer embedded between Al- 
GaAs barriers [3(| • This part of the paper is intended to 
demonstrate the potential of our theory within another 
application and to stimulate further experiments for this 
material system. In a good approximation, the carriers 
in the GaAs layer are confined by infinitely high potential 
barriers in the z-direction . The solid line in Fig. [5^ 
shows the linear optical transmission T through a GaAs 
layer of 5 exciton Bohr radii thickness. The energy fko is 
given relative to the bulk band-gap energy Eg and nor- 
malized to the corresponding exciton binding energy Eg . 
For the chosen layer thickness, the confinement of carri- 
ers in the z-direction yields three polariton resonances 
of the Is exciton in the displayed part of the spectrum 
which are labeled by consecutive numbers. 

For the pump-probe excitation, two 120 fs laser pulses 
(from slightly different directions) without time-delay are 
considered. The dashed line in Fig. [5^ corresponds to 
the pulse spectrum. Figure |2b shows changes in the 
probe-pulse transmission that are induced by the pump 
pulse for opposite circular (solid line) and co-circular 
(dashed-dotted line) polarization of pump and probe 
pulses. The probe pulse enters the excitonic polariza- 
tion in linear order only. For the pump pulse a Rabi en- 
ergy d c h|E pump | = 0.01 Eg is used which ensures a consis- 
tent description within our theory. The transmission 
changes around the higher (Is, 2 and Is, 3) polariton re- 
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FIG. 2: (a) Solid line: Calculated linear transmission spec- 
trum T (depicted as 1 — T) for a GaAs layer with 5 exciton 
Bohr radii thickness. Dashed line: Spectral shape of a 120 fs 
laser pulse, (b) Differential probe transmission AT for oppo- 
site circular (solid line), and co-circular (dashed-dotted line) 
polarization of pump and probe pulses, (c) Differential probe 
transmission AT for opposite circular polarization including 
all Coulomb terms (solid line) and without Coulomb interac- 
tion of different polaritons (dashed line) . 



sonances in Fig-Et* are similar to those around the lowest 
one (ls,l) but with a decreased amplitude. A similar de- 
pendence on the light polarization has been reported for 
the differential probe absorption around the Is exciton 
resonance in a QW system |22( . For opposite-circular po- 
larization, the pump-induced changes in the probe trans- 
mission are exclusively determined by biexcitonic corre- 
lations no mean- field effects contribute. This con- 
figuration is chosen for the analysis of the Coulomb in- 
teraction between polaritons in states with different spa- 
tial distribution. The dashed line in Fig. shows the 
result where Coulomb interaction that couples different 
eigenf unctions in the two-exciton product basis of Eq. (@J 
is artificially switched off. We encounter only a slight 
change of the nonlinearities around the lowest polariton 
resonance (ls,l) whereas for higher peaks (Is, 2 and Is, 3) 
the influence of the pump pulse almost vanishes. There- 
fore, Coulomb interaction between exciton states with 
different spatial distribution (corresponding to different 
polariton resonances) turns out to be the main source for 
transmission changes of higher polariton states. 

In summary, a theory for the nonlinear polariton dy- 
namics in semiconductor heterostructures has been pre- 
sented which consistently includes (i) propagation effects 
with microscopic boundary conditions for the induced 
material polarization as well as for the optical fields 



and (ii) excitonic and biexcitonic coherent nonlineari- 
ties previously studied only in QWs or one-dimensional 
model systems. The application to nonlinear transmis- 
sion spectra for a ZnSe/ZnSSe heterostructure shows ex- 
cellent agreement between theory and experiment. The 
influence of light-polarization dependent excitonic and 
biexcitonic nonlinearities has also been demonstrated 
for a pump-probe excitation of a GaAs heterostructure. 
Coulomb interaction of polariton states with different 
spatial distribution turns out to strongly influence the 
nonlinear optical transmission. 
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